Ballistic (precessional) contribution to the conventional magnetic switching 
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We consider a magnetic moment with an easy axis anisotropy energy, switched by an external 
field applied along this axis. Additional small, time-independent bias field is applied perpendicular 
to the axis. It is found that the magnet's switching time is a non-monotonic function of the rate at 
which the field is swept from "up" to "down". Switching time exhibits a minimum at a particular 
optimal sweep time. This unusual behavior is explained by the admixture of a ballistic (precessional) 
rotation of the moment caused by the perpendicular bias field in the presence of a variable switching 
field. We derive analytic expressions for the optimal switching time, and for the entire dependence 
of the switching time on the field sweep time. The existence of the optimal field sweep time has 
important implications for the optimization of magnetic memory devices. 



In conventional magnetic switching by an externally 
applied magnetic field the moment performs many rev- 
olutions before switching to the opposite direction. 
While being much slower than the ballistic (precessional) 
switching^— which is tested experimentally but not yet 
realized in applications, conventional switching is used 
for magnetic recording in hard disk drives and other de- 
vices. The speed at which the moments of magnetic bits 
can be switched between the two easy directions has obvi- 
ous implications for the technology performance, setting 
the limit for the information writing rate. Here we study 
the dependence of the conventional magnetic switching 
time r m on the reversal time of the writing head field. If 
the field is swept from "up" to "down" in a time r^, the 
switching time will be a function T rn (Th). It would seem 
natural to assume that r m decreases with decreasing t% 
and the fastest switching is realized by an instantaneous 
flip of the field with = 0. However, it was found nu- 
merically by one of the authors^ that in the presence of 
a small perpendicular bias field the function T m (rh) is 
not monotonic and has a minimum at an optimal sweep 
time t?. Decreasing Th below the optimal value would 
be counterproductive in terms of the technology perfor- 
mance. In this paper we provide analytic approximations 
for the function r m (r/j) and the optimal field sweep time 
Ty t . We find that the nonmonotonic behavior of r m {jh) 
is a result of the admixture of a "ballistic" (or "pre- 
cessional") switching induced by the perpendicular bias 
field. Ballistic contribution is normally quenched by the 
anisotropy, but here it is restored by the time-dependence 
of the field during the rise time of the applied step. The 
effect considered here is different from the decrease of 
the switching time^ predicted for switching below the 
Stoner-Wholfarth limit The latter consists of the r m 
dependence on the amplitude of the field step with an 
instantaneous rise time, while in our case the finite rise 
time is essential. Conventional switching considered here 
is achieved by a field step and does not require precisely 
timed pulses of finite duration needed for a truly ballistic 
switching. 1-3 Also, the required field magnitude is much 
smaller than in the ballistic case. 



Magnetic particles of nanometer size are single- 
domain 1 '^ and can be described by the moment M = 
M n, where n is a unit vector. We consider a par- 
ticle with an easy axis z and anisotropy energy E a — 
-(l/2)Kn 2 z . The switching field H = H(t)z is directed 
along the easy axis. For large enough field magnitudes, 
|.ff | > K/Mo, only one equilibrium direction of M is sta- 
ble. A field applied exactly along the axis leads to a 
magnetization switch only when some fluctuations of M 
are present. Following Ref. |4j we introduce a small con- 
stant perpendicular bias field = H±x, H±_ -C K/Mo 
to mimic the required fluctuations. 

The dynamics of the moment are governed by the 
Landau-Lifshitz-Gilbert (LLG) equation 
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where 7 is the gyromagnetic ratio, E — E a — (H+H^)-M 
is the total magnetic energy, and a <tc 1 is the Gilbert 
damping constant. The field sweep is assumed to be lin- 
ear in time and given by the expressions H(t) = +Hq 
for t < 0, H(t) = H (l - 2t/r h ) for < t < r h , and 
H{t) = —Hq for t > 77,. As the field is swept from pos- 
itive to negative values, the up-equilibrium disappears 
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FIG. 1: Switching time r m as a function of field sweep time 77,. 
Here a = 0.01, ho = 3.5 ljo, h± = 0.001 cjo- The solid line is 
given by the analytical expression {TJ. 
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FIG. 2: Dependencies T m {jh) calculated for ho = 2.2 Wo, 
ftx = 0.001 cjq, and variable a indicated on each panel. As a 
increases, the theoretical fit gets poorer due to the violation 
of the strong inequality (j4j). 



and the magnetic moment starts to move towards the 
down-equilibrium along a spiral trajectory. The final ap- 
proach to the down-equilibrium in exponential. To define 
a finite switching time we have to introduce a provisional 
cut-off angle 9 SW and calculate the time it takes to reach 
it. The extra time needed to cover the remaining dis- 
tance does not depend on 77, because that part of the 
motion happens at a constant field H = —Hq. In accord 
with Refs. |j-0 we choose 9 SW = ir/2. 

The LLG equation can be easily solved numerically and 
the switching time dependence T m {jh) can be obtained. 
Fig.[T]shows the results of such modeling for a particular 
parameter set. The minimum of r m is clearly observed. 

We have derived an approximate analytic expression 
for the switching time. Denoting ho — jHq, h± = 7-Hj_, 
and wq = jK/Mq we find 
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Formula ([1} is the first main result of our paper. As one 
can see in Fig.[TJ it reproduces the function T m {jh) quite 
well. 

Approximation (JT|) requires small h± and a. For a 
given h± it is valid in the interval of field sweep times 
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where is a solution of 
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The optimal sweep time is determined from 
dT m /dTh = 0. We get an expression 
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This formula is our second main result. Note that t! 
is independent of the bias field. The minimal switching 
time T m (r^) itself depends on h±, which is quite natural 
since the initial deviation from the easy axis is controlled 
by h±. We have also calculated the switching time drop 
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between the instant and the optimal field sweeps. The 
drop is independent of the bias field as long as approxi- 
mation ([lj is valid. 

Formula © is valid when falls into the interval © . 
Our calculations show that this is guaranteed for 
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These inequalities place a more stringent constraint on 
the Gilbert damping than the simple a <C 1. 

Fig.[l]compares numerically calculated switching times 
with our analytic formula (H). When inequalities (j3| are 
well satisfied, the quality of approximation is very good. 
As one approaches the limits of the approximation's va- 
lidity by, e.g., increasing a, the errors grow larger. 

Fig. |3] shows the dependence of the optimal field sweep 
time on the system parameters. The correspondence 
with formula <j3j> is generally good, although some visible 
deviations exist. The accuracy of the determination of 
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FIG. 3: Numeric (points) and approximate analytical (solid 
lines) dependencies of the optimal field sweep time on the 
system parameters: (a) fixed a and h±, (b) fixed ho/u>o and 
h±, (c) fixed fto/wo and a. When not varied, the parameter 
values are ho = 2.2wo, a = 0.01, h± = 0.005. 
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FIG. 4: Average "ballistic" contribution of the bias field. 
Vector M(£) orbits around a parallel circle C on a sphere 
|M| = Mo. The torque due to Hi pushes M along the 
meridians of the sphere. In constant switching field H the 
torque contributions from the diametrically opposed elements 
dh and dh cancel each other. For variable U(t) such a cance- 
lation does not happen because M spends unequal amounts 
of time on dli and dfa- 

t£ is lowered by a flat shape of the r m (r^) minimum. 
The shallow minimum, however, also lowers the practical 
importance of precise determination of . 

In general, the analytic expression can approximate the 
T m (Th) dependence up to a 10% accuracy in a surprisingly 
wide range of parameters. Such accuracy is certainly 
sufficient for the estimates related to the device design. 

We now discuss the physical reason for the minimum 
of the function T m (T/j). The bias field has two roles in the 
switching process. First, it provides the initial deviation 
from the easy axis. Second, it alters the equations of mo- 
tion for M(i). The derivation of Eqs. (JlJ and which 
will be detailed in our forthcoming paper, shows that the 
second contribution is dominant. Recall now that in the 
absence of anisotropy and other fields the torque H± x M 
due to the bias field would rotate vector M from +z to 
— z along a meridian of the sphere |M| = M$ (dashed line 
in Fig.[4j), in a ballistic (precessional) fashion. In our case 



a weak bias field is applied on top of the strong uniax- 
ial anisotropy and switching field, which together induce 
a fast orbital motion of vector M(i) along the parallel 
circles (line C in Fig. @|. The bias field still attempts 
to move M along the meridians, but now its action has 
to be averaged over the orbital period. As illustrated in 
Fig. |4j in constant fields H = ±Hqz averaging gives zero 
due to the cancelation of the contributions from the di- 
ametrically opposed infinitesimal intervals dl\ and dli of 
equal lengths. This way ballistic contribution of the bias 
field is quenched. However, the contribution of Hj_ does 
not average to zero for a variable switching field H(t). In 
this case the velocity of M changes along the orbit, the 
times spent in the intervals dl\ and dli are different, and 
the contributions of the two do not cancel each other. We 
conclude that in the presence of a time dependent exter- 
nal field H(i) ballistic contribution of the perpendicular 
bias field is recovered. Moreover, this contribution helps 
to move vector M from +z to —z and is thus respon- 
sible for the initial decrease of r m . As the sweep time 
grows larger, the change of the orbital velocity during the 
precession period decreases and the ballistic contribution 
averages out progressively better. The helping effect of 
ballistic switching is lost and T m starts to increase as it 
normally would. 

Ballistic contribution to switching can be also viewed 
as a phenomenon complimentary to the magnetic reso- 
nance and rf-assisted switching,— ~— where H is constant 
but Hj_(£) is time-dependent. There the average contri- 
bution of the bias field on an orbit does not vanish due to 
the time dependence of H±. The non- vanishing contri- 
bution, regardless of its origin, assists the switching and 
makes it faster. 

Our analytical results provide a convenient approxi- 
mation for the optimal field sweep time, an important 
parameter in the device design. They can be used as a 
starting point for the investigations of the switching time 
in granular media, where each grain can be modeled by a 
single moment and bias fields are produced by the other 
grains or by the spread of grain orientations. 
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